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^ '. ABSTRACT 

in . 
o ■ 

\ H1426+428 is one of blazars that are observed by 7-rays in the TeV region. 

P-f Because TeV 7-rays from distant sources are subject to attenuation by the ex- 

tragalactic background light (EBL) via electron-positron pair production, the 
intrinsic spectrum of the TeV 7-rays should be inferred by using the models of 
radiation processes and EBL spectrum. We set constraints on the physical con- 
.^J I dition of H1426+428 with the synchrotron-self-Compton model applying several 

EBL models. We find that the emission region of H1426+428 is moving toward 
us with the bulk Lorentz factor of ~ 20 and that its magnetic field strength is 
~ 0.1 G. These properties are similar to other TeV blazars such as Mrk 421 and 
Mrk 501. However, the ratio of the energy density of nonthermal electrons to 
that of the magnetic fields is about 190 and fairly larger than those of Mrk 421 
and Mrk 501, which are about 5 - 20. It is also found that the intensity of EBL in 
the middle and near infrared wavelengths should be low, i.e., the intensity at 10 
fim is about 1 nW m~ 2 sr _1 to account for the observed TeV 7-ray flux. Because 
the spectral data of H1426+428 in X-rays and 7-rays used in our analysis were 
not obtained simultaneously, further observations of TeV blazars are necessary 
to make the constraints on EBL more stringent. 



Subject headings: BL Lacertae objects: individual (H1426+428) - diffuse radi- 
ation - galaxies: active - radiation mechanisms: nonthermal - infrared: general 
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1. Introduction 

Recent observations of very high energy 7-rays revealed that some blazars emit 7-rays in 
the TeV region (e.g., Aharonian et al. 2005b, and references therein). The spectral energy 
distribution (SED) of these blazars has two peaks; one in the X-ray region and the other in 
the very high energy 7-ray region (GeV - TeV). The emission mechanisms of TeV blazars 
are thought to be synchrotron radiation and inverse Compton scattering by nonthermal 
particles in relativistically moving jets (e.g., Maraschi et al. 1992; Marscher & Travis 1996; 
Inoue & Takahara 1996; Li & Kusunose 2000; Kino et al. 2002). Leptonic models that 
assume electrons and positrons as emitting particles are most widely accepted. Although 
hadronic models are viable, which postulate emission from mesons and leptons produced 
by the cascade initiated by proton-photon or proton-proton collisions (e.g., Mannheim & 
Biermann 1992; Mannheim 1993; Miicke et al. 2003), we do not consider these models in 
this paper. 

It has been known that TeV 7-rays emitted by blazars are attenuated by absorption 
due to the extragalactic background light (EBL) en route to the Earth, because e + e~ pair 
production occurs via 7 + 7 — * e + + e - . The importance of the absorption of high energy 
7-rays due to cosmic background radiation was noticed quite early (Nikishov 1962; Gould 
& Schreder 1966; Jelly 1966). The effects of the attenuation on the SED of blazars were 
recognized by Stecker et al. (1992) following EGRET observation of 3C 279 (Hartman et 
al. 1992). Since then this effect has been extensively studied by many authors (e.g., Dwek 
& Slavin 1994; Stecker & Salamon 1996; Konopelko et al. 2003; Kneiske et al. 2004; 
Dwek & Krennrich 2005). The spectrum of EBL (A ~ 0.1 to 1000 //m) has been observed by 
various instruments (Hauser & Dwek 2001; Kashlinsky 2005, for review). It is known that 
there are two peaks in the SED of EBL; one at A ~ 1 jum is associated with star light and 
the other at A ~ 100 /im is related with dust emission. The wavelength region of ~ 5 to 60 
/im is difficult to observe because of interplanetary dust emission, while the attenuation of 
TeV 7-rays is most effective in this wavelength region of EBL. When a 7-ray with energy 
collides with an EBL photon with energy £ebl, pair production occurs if E^Eebl > ijn e c 2 ) 2 
is satisfied, where m e and c are the rest mass of an electron and the light speed, respectively. 
If E y = 1 TeV, a photon with Aebl = he/ £ebl < 1 A*ni produces a e + e~~ pair, where h is the 
Planck constant. 

Because both the intrinsic SED of TeV blazars and the EBL spectrum are not well 
known, various attempts have been made to construct self-consistent models of TeV blazars 
and the EBL spectrum. (Differences in various EBL models are discussed in detail by 
Hauser & Dwek (2001).) For example, De Jager & Stecker (2002) derived the intrinsic 
7-ray spectrum of Mrk 501 during 1997 high state, using the models of EBL by Malkan & 
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Stecker (2001) (backward evolution model). Recently Dwek & Krennrich (2005) employed 
a different way of constructing the models of the EBL spectrum. They adopted various 
observed EBL data and fitted the data with polynomials; 12 different EBL template spectra 
with different intensities in UV - near IR, mid IR, and far IR spectra. Using those spectra, 
they set constraints on the intrinsic spectra of Mrk 421, Mrk 501, and H1426+428 and the 
acceptable range of the EBL spectrum was estimated. 

Most models of the SED of TeV blazars and the EBL spectrum have not taken into 
account the electron spectrum in the emission region self-consistently. Many authors assumed 
a given form of the energy spectrum of electrons or do not consider the emission mechanisms 
at all. Konopelko et al. (2003), on the other hand, solved kinetic equations of electrons 
and photons simultaneously to calculate the emission spectrum of Mrk 501 and Mrk 421. 
They calculated the optical depth for TeV 7-rays according to the EBL models by Malkan 
& Stecker (2001) and De Jager & Stecker (2002). They obtained a fairly good fit to the 
observed spectrum of the blazars with a rather large value of the Doppler factor T> of the 
jet, i.e., V ~ 50. 

In this paper we make models of H1426+428 (at redshift z = 0.129) and estimate 
the intensity of EBL, solving the kinetic equations of electrons and photons to obtain the 
intrinsic SED of blazars, so that the observed SED is consistent with the electron injection 
and cooling processes. X-ray observation of H1426+428 was carried out by Costamante et 
al. (2001) with BeppoSAX and X-ray spectrum was found to extend up to 100 keV. Falcone 
et al. (2004) observed H1426+428 with RXTE and found that the peak energy in SED 
was sometimes in excess of ~ 100 keV and at other times in the 2.4 - 24 keV region. TeV 
7-rays have been observed by several groups (Horan et al. 2002; Aharonian et al. 2002; 
Djannati-Atai et al. 2002; Petry et al. 2002; Aharonian et al. 2003). Because the redshift 
of H1426+428 is 0.129, the effects of the attenuation of TeV 7-rays should be very large and 
the object can be a good test against various EBL models. 

Costamante et al. (2003) inferred the intrinsic spectrum of H1426+428 using the model 
spectrum of EBL given by Primack et al. (2001), who employed the semi-analytical models 
of galaxy formation, as well as a EBL model spectrum designed to match the high fluxes 
below 2 /im (Aharonian et al. 2002). Their model spectrum of H1426+428 has a large 
value of the ratio of Compton luminosity to synchrotron luminosity, Lq/ L s > 10. According 
to Costamante et al. (2003), external seed radiation fields for inverse Compton scattering 
are necessary to account for Lq/L s > 10 and to avoid the Klein-Nishina effects in inverse 
Compton scattering. Aharonian et al. (2003) obtained a rising intrinsic spectrum in the 
TeV regime, using the same EBL spectra as Costamante et al. (2003) and Malkan & Stecker 
(2001) adopted. 
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Since the estimate of the intrinsic SED of TeV blazars largely depends on the EBL 
spectrum and the uncertainty of the SED of EBL is quite large, it is worth testing various 
possibilities of the EBL spectrum. In this paper we use the templates of EBL derived by 
Dwek & Krennrich (2005), because their model spectra cover a wider range in intensity 
than those used by Konopelko et al. (2003). We present the method of calculations of the 
energy spectra of photons and electrons in the jet and the spectrum of EBL in §2. Numerical 
results for energy spectra of photons and electrons in the emission region of H1426+428 are 
shown in §3. Finally summary of our results and discussion are given in §4. 



2. Model 

2.1. Kinetic equations 

The emission mechanisms of H1426+428 are assumed to be synchrotron radiation and 
synchrotron-self-Compton (SSC) scattering. We calculate the SED of H1426+428 before 
the attenuation of 7-rays by simultaneously solving the kinetic equations of photons and 
electrons in the emission region. We use the same numerical scheme as used in Kusunose 
et al. (2003) except that external soft photons are not included in this work. We assume 
that the emission region is a spherical blob with radius R moving with a relativistic speed 
with Lorentz factor V. It is also assumed that T> — V. Nonthermal electrons are injected 
into the emission region at injection rate q e (l) per unit 7 and unit volume in the comoving 
frame of the blob, where 7 is the electron Lorentz factor. The kinetic equation of nonthermal 
electrons is given by 

^M + ^M = ^ [ftyn+ ^ oK(T)] + , e(T ), (1) 

where t eiesc is the electron escape time from the emission region, and 7 syn and 7ssc are the 
cooling rates by synchrotron radiation and SSC scattering, respectively. Here 7ssc depends 
on the photon spectrum in the blob. The injection rate of nonthermal electrons is given by 

Qe(l) = %1~ P exp(-7/7 max ) , 7 > 7 min , (2) 

where p, 7 m i n , and 7 max are parameters, which are determined by fitting the emission spec- 
trum with observed one. 

The kinetic equation of photon number spectrum n 7 (e), where e is the photon energy, 
is given by 

^M + ^lM = n c (e) + n syn {e) - h ahs (e) , (3) 

OZ £7, esc 
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where i 7)esc = R/c is the photon escape time from the emission region. The photon produc- 
tion rate per unit volume in the blob by synchrotron radiation is denoted by n syn (e). The 
self-absorption by synchrotron radiation and the absorption by e + e~ pair production inside 
the blob are included in n a bs(^)- The production and loss of photons by Compton scattering 
is given by 



h c (e) = -n J (e) / ^7^(7) R c (e, 7) + / / de' d-f P(e;e', 7) Rc{e' ,l) n 7 (e')n e (7) , (4) 



where Rc(s, 7) is the angle averaged scattering rate and P(e;e', 7) is the probability that 
photon energy e' is changed to e by a single scattering with an electron with 7 (Coppi & 
Blandford 1990). Here the exact Klein-Nishina cross section is used to calculate Rq and 
P. We solve equations (1) and (3) simultaneously to obtain the spectra of electrons and 
photons self-consistently. The emission spectrum is calculated from the term of photon 
escape in equation (3) and the observed spectrum is affected by the Doppler and beaming 
effects because of the relativistic motion of the emission region. We include these effects to 
compare the model SED with observations. 



The interaction between high energy 7-rays and EBL via e + e _ pair production attenu- 
ates TeV 7-rays from blazars. The optical depth of the attenuation r 77 is a function of the 
7-ray energy and the EBL spectrum. For the calculation of r 77 we assume a flat ACDM 
cosmology with H = 71 km s" 1 Mpc -1 , Vt m = 0.27, and fi A = 0.73 (Bennett et al. 2003), 
where H is the Hubble constant and Vt m and are the density parameters of matter and 
dark energy, respectively. The EBL number density at z is calculated as 



where no(u ) is the EBL number spectrum measured at z = and u = v(\ + z)~ l . Equation 
(5) assumes that the phase space density of EBL is conserved. The spectrum of EBL has been 
measured by various groups. However, there is a large uncertainty in the observed spectra. 
To take into account various possibilities, we follow the recently adopted method by Dwek 
& Krennrich (2005). They generated 12 EBL templates based on numerous observational 
data. The templates are labeled as XYZ such as LLL, MLH, and LHH, etc. Namely, X = 
L, M, or H representing the energy flux of the stellar component of EBL (L = low, M = 
medium, and H = high). Also, Y = L or H, and Z = L or H, where Y denotes the flux level 
at 15 (Ltm band and Z does at far-IR band. 




2.2. Optical Depth by EBL 



^ebl(^, z) = (1 + z) 2 n (u ) , 



(5) 
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The observed SED is calculated by multiplying the model SED by exp(— r 77 ). Because of 
this exponential factor of the attenuation, the differences in r 77 are significant in construction 
of the SED models of blazars. 



3. Numerical Results 

In Figure 1, the 7-ray data obtained by the Whipple collaboration (Horan et al. 2002) 
and the HEGRA collaboration in the term 1999 - 2000 (Aharonian et al. 2003) are shown. 
These are calculated by multiplying the observed data by exp(r 77 ) assuming various models 
of EBL. It is found that the attenuation free SED does not depend much on Z of EBL, i.e., 
SED is almost the same for Z = L and Z = H. To distinguish the dependence on Z, 7-ray 
sources with higher energies are needed. When EBL models with Y = L are applied, SED 
is flat or becomes downward as the photon energy increases. For EBL models with Y = H, 
the peak of vF v before the attenuation may appear at photon energies greater than 50 TeV, 
which is hard to explain by SSC models; this will be discussed in §4. Because we search for 
the values of parameters with SSC models, EBL models with Y = H are excluded. 

Our model SEDs of H1426+428 are shown in Figure 1 by solid and dashed lines. The pa- 
rameters to specify the spectra of electrons and photons of the emission region of H1426+428 
are T, R, magnetic field B, p, 7 m i n , 7 max , %, an d t e)CSC . The model spectrum of EBL that 
can be consistent with our SED models are LLL, LLH, MLH, and MLL. With these EBL 
templates, it is not necessary to invoke the external soft photons as seed photons of inverse 
Compton scattering, because of the smaller value of the optical depth for 7-rays in 1 - 20 
TeV region. We first adopt X = L as a template for which the model of H1426+428 is 
obtained. In Figure 2, the 7-ray data together with the data in other wave bands are shown, 
where HEGRA data in the term 2002 are also shown for reference. The X-ray data obtained 
by BeppoSAX are read from Figure 3 in Costamante et al. (2003) and only three points 
are plotted in this figure. It is to be noted that the data in multiwavelength bands were 
not obtained simultaneously. We find that the following values of the parameters for the 
model SED before the attenuation fit the observed data: B = 0.1G, r = 20, R = 10 16 cm, 
p = 1.9, 7 m i n = 2, 7 max = 10 7 , the electron injection rate of 10~ 3 electrons cm~ 3 s _1 , and 
^e,esc = 4i?/c. We call this model "base model," hereafter. (The solid line in Figure 1 is the 
SED calculated with these parameters.) In Figure 2, the attenuation free SED is shown by 
a dashed line. The model SEDs in the observer frame are calculated by multiplying the base 
model by exp(— r 77 ) with the EBL models, LHH and LLL. The EBL models LHH and LHL 
result in almost the same SEDs, and LLL and LLH also result in almost identical SEDs. 
In these calculations, we included the pair production that occurs in the emission region, 
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but this effect is found to be negligible. The peak value in the v-vF v representation of the 
SSC component is given at hv ~ 8.5 x 10 10 eV and ~ 6.5 x 10 10 eV for before and after the 
attenuation by EBL, respectively. These values are much smaller than those obtained by 
Costamante et al. (2003). Also these are smaller than the peak energy obtained for Mrk 
421 and Mrk 501 by Dwek & Krennrich (2005). 

For the base model shown in Figure 2, the energy contents of the emission region are as 
follows; magnetic energy density M mag = 3.98 x 10~ 4 ergs cm" 3 , nonthermal electron energy 
density u e = 7.51 x 10 -2 ergs cm -3 , and photon energy density w p h = 8.66 x 10~ 2 ergs cm~ 3 . 
Thus « e /« mag R3 188 and this value is larger than u e /u ma , g ~ 10 for other TeV blazars (Kino 
et al. 2002). This large value of M e /« mag is a result of a large energy density of very high 
energy electrons, which is necessary to emit a large flux of TeV 7-rays. The high energy 
tail of the 7-ray spectrum is highly affected by the Klein-Nishina effects in inverse Compton 
emission and the attenuation by EBL. If the flux of EBL is larger than LLL or LLH model, 
much larger values of u e and V are needed to account for the observed flux of TeV 7-rays. 

The energy spectrum of nonthermal electrons is shown in Figure 3 for the base model. 
The injected electron spectrum has an exponential cutoff at 7 = 10 7 . The steady state 
energy spectrum of electrons has a peak at 7 ~ 10 4 that is caused by the balance between 
radiative cooling and escape of electrons from the emission region. The observed character- 
istic synchrotron frequency due to electrons with 7 = 10 4 is about 8 x 10 14 Hz, if V = 20 
and B = 0.1 G. The value of 7 m j n is 2 and the cooling time for electrons with this value of 7 
is so long that there is no flow in the energy space of electrons toward lower energies below 
7 = 2. 

In Figure 4, a model SED with X=M, i.e., MLL, is shown. To fit the data with MLL, a 
very large SSC component is necessary. For this, B = 0.02 G and the injection rate of 5 x 10~ 3 
electrons cm~ 3 s _1 are assumed, and other parameter values are the same as those of the base 
model. (The dashed line in Figure 1 is the SED calculated with these parameters.) Because 
the SSC component highly dominates over the synchrotron component, an extremely large 
value of u e /u mag is realized, i.e., ~ 2.87 x 10 4 . 

We compare the models of the emission spectrum of H1426+428 with different values of 
parameters in Figure 5. Note that in this figure observed 7-rays are corrected by multiplying 
exp(r 77 ) to reproduce the SED before the attenuation, assuming LLL model for EBL. When 
the magnetic fields are weaker but the other parameter values are the same as for the base 
model in Figure 2, the synchrotron cooling rate decreases and the number of high energy 
electrons increases. As a result the SSC component has a larger luminosity. An example 
with B = 0.05 G is shown by a dotted line. When the value of 7 m i n is larger, the synchrotron 
luminosity in radio decreases and the inverse Compton luminosity in the MeV region also 
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decreases. In Figure 5, 7 min = 10 3 is shown by a dash-dotted line. This SED is calculated 
with the value of q adjusted so that the peak luminosity in X-rays is almost the same as 
that of the base model. The difference in t e esc does not change much the spectra of X-ray 
and TeV 7-ray, although the larger value of t eiesc results in a larger flux in MeV 7-rays. The 
important parameters to determine the TeV flux are 7 max and T. When the value of 7 max is 
larger, the peak energy of synchrotron emission becomes larger, but the peak energy of the 
SSC component does not change much because of the Klein-Nishina effects. On the other 
hand, the value of T does not change the spectral shape. It appears that the base model 
does not yield a good fit in X-ray regime, and the allowable parameter ranges are still large. 
Since the data in X-rays and 7-rays were not obtained simultaneously stronger constraints 
are difficult to set. 



4. Summary and Discussion 

We obtained the self-consistent models of the emission region of H1426+428, taking 
account of the interaction of TeV photons and EBL photons via 7 + 7 — > e + + e~. In 
the calculations of the emission spectrum of H1426+428, we solved the kinetic equations 
of electrons and photons in a relativistically moving blob. To obtain the optical depth for 
the attenuation (7 + 7 — > e + + e~) we used 12 templates of EBL constructed by Dwek & 
Krennrich (2005). 

We found the parameter sets of H1426+428 for the models with EBL models with X = 
L or M and Y = L. Here X and Y represent the energy fluxes of the stellar component and 
15 /im band, respectively. Other EBL models require stronger emission of TeV 7-rays and 
we did not find better models of H1426+428. With Y = L, we found that the bulk Lorentz 
factor of H1426+428 is about 20 and that the magnetic field is ~ 0.1 G. These are similar 
values to those of other TeV blazars such as Mrk 421 (Kino et al. 2002). The energy density 
of the nonthermal electrons dominates that of magnetic fields. This is the same properties 
as have been found for other TeV blazars such as Mrk 421 and Mrk 501 (Kino et al. 2002), 
but the value of w e /w mag is about 190 and this is much larger than the values of about 10 
for Mrk 421 and Mrk 501. The injected electron number spectrum follows a power law with 
a power-law index of 1.9 and has an exponential cutoff at 7 max = 10 7 . This value of 7 max is 
larger than that for Mrk 421 and Mrk 501 (Kino et al. 2002) and very efficient acceleration 
mechanisms is suggested in this source. The electron energy spectrum in the steady state in 
the emission region has a break at 7 ~ 10 4 . Very efficient cooling by synchrotron radiation 
and inverse Compton scattering produces this break and a large fraction of particle energy 
is lost as radiation. In many models to explain the high energy 7-rays from blazars, the 



- 9- 



electron spectrum in the emission region is just assumed to obey a power law or a broken- 
power law. However, it is clear from our results that solving the kinetic equation of particles 
is necessary to obtain the self-consistent particle spectrum in particular when Klein-Nishina 
effects are significant. 

The peak value of vF v in the 7-ray region is attained at his ~ 8.5 x 10 10 eV before the 
attenuation according to our base model. This value is significantly smaller than the peak 
energy estimated by Costamante et al. (2003) for H1426+428 (larger than 10 TeV) and 
the values estimated for Mrk 421 and Mrk 501 (~ 1-3 TeV) by Dwek & Krennrich (2005). 
Costamante et al. (2003) concluded that inverse Compton scattering of external soft photons 
produces TeV 7-rays, because SSC models suffer from the Klein-Nishina effects. However, 
our results show that there is no need for the external soft photon sources. This discrepancy 
originates from the use of different models of EBL. Costamante et al. (2003) used the 
models by Aharonian et al. (2002) and Primack et al. (2001) and these models yield 
higher intensity of EBL than LLL and LLH. Then the similar arguments might be raised, 
if EBL models such as HHH, HHL, MHH, and MHL are applied. That is, for these EBL 
models, TeV emission appears to have a peak above 50 TeV in SED before the attenuation 
and inverse Compton scattering of external soft photons might be necessary to account for 
the 7-rays spectrum. If the X-ray emission at ~ 10 19 Hz (a peak frequency in SED) is by 
synchrotron emission, radiating electrons have Lorentz factor of ~ 10 6 , where B = 0.1 G 
and T = 20 are assumed. To produce 7-rays with energy of 50 TeV via inverse Compton 
scattering by the same electrons in the Thomson regime, the energy of soft photons in the 
rest frame of external soft photon sources is in the IR region with the frequency of about 
10 13 Hz and the broad line regions can be a photon source. Because of the large cooling 
rate of Compton scattering of external soft photons, the values of B, 7 max) and go will be 
modified. It is, however, beyond the scope of this paper to find another model with Compton 
scattering of external soft photons. 

Dwek k Krennrich (2005) analyzed the data of Mrk 421, Mrk 501, and H1426+428, 
and they rejected some EBL spectra which cause "unphysical intrinsic spectra characterized 
by an exponential rise at high TeV energies." The EBL models, LLL, LLH, MLL and MLH, 
which are found to be good candidates for the EBL spectrum in our model, are also allowed 
models in Dwek & Krennrich (2005). Since they only tested the behavior of high energy 7- 
rays, their constraints may not be restrictive. In contrast, we searched for a set of parameter 
values that realizes multiwavelength properties of H1426+428, based on the homogeneous 
SSC model. 

In this paper we assumed that the phase space density of EBL is constant [equation (5)]. 
It is possible that the EBL spectrum may not change with redshift so that ^ebl(^5 z)du = 
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no{vo)dvo, because stellar population and dust emission might be almost homogeneous for 
small values of z. We tested this model of EBL and found the value of r 77 of LLL at E — 1 
TeV decreases about 18 %. Then above mentioned results with the constant phase space 
density of EBL are not changed much. 

Because the multiwavelength data of H1426+428 have not been obtained simultane- 
ously, the constraints on the physical condition of H1426+428 and the EBL spectrum are 
not stringent yet. Also, because there are several parameters to determine the emission spec- 
trum from the jet, the allowed ranges of the parameter values are not strictly determined. 
According to our numerical calculations, the magnetic fields, for example, might be in the 
range between 0.05 and 0.2 G, and the bulk Lorentz factor might be in the range between 
15 and 25. However, the properties of H1426+428 found in this work will not be changed 
dramatically as long as the SSC model is applied and the EBL flux at around 10 /im should 
be as low as the model values such as LLL, LLH, MLL, and MLH. 

Recent observations of PKS 2155-304 by H.E.S.S. (High Energy Stereoscopic System) 
showed that it emits 7-rays above 160 GeV and possibly above a few TeV (Aharonian et al. 
2005a). Because this object is located at z — 0.116, the absorption of TeV 7-rays by EBL 
is probable and further constraints on the EBL spectrum will be obtained. More recently, 
Aharonian et al. (2005c) reported the detection of TeV 7-rays by H.E.S.S. from further 
distant blazars, i.e., H2356-309 (z = 0.165) and 1ES 1101-232 (z = 0.186). Assuming 
that the source spectra of 7-rays obey a power law, dN/dE oc E~ a , where dN/dE is the 
differential number spectrum of 7-rays and a is the power-law index, they set upper limits 
on the EBL spectrum in the wavelength band ~ 0.1 - 10 /im. The EBL models of LLL and 
LLH that are favored in this work are consistent with the EBL models of P0.45 and P0.55 
of Aharonian et al. (2005c) in A = 1 - 10 /xm. 

Dwek et al. (2005) recently examined the extragalactic nature of the excess near- 
infrared background light (NIRBL) (e.g., Matsumoto et al. 2005), using the analysis of the 
7-ray spectra of blazars such as H1426+428 and PKS 2155-304. This excess NIRBL is still 
controversial and they argued that the excess NIRBL may not be related with the spectral 
imprint of the first generation of stellar objects on the EBL. 

In this paper we assumed the steady state of the emission region, but the time variation 
of X-rays and TeV 7-rays of blazars is common. The time evolution of the emission spec- 
trum might be used to constrain the emission mechanisms of blazars including external soft 
photons and hadronic models. Then the EBL spectrum might be constrained further. Since 
our numerical scheme is capable of time-dependent calculations, we may apply our model to 
time variable blazar emission in future work. 
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Fig. 1. — TeV spectrum [data from Whipple (Horan et al. 2002)] and HEGRA [data from 
(Aharonian et al. 2003)], corrected by 12 different EBL models. The solid lines are our 
model SED before the attenuation (base model), which is shown by a dashed line in Fig. 2. 
The dashed lines are also our model SED before the attenuation shown by a dashed line in 
Fig. 4. 
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Fig. 2. — SED of H1426+428 and model spectra with different EBL models. Data obtained 
in 2002 by HEGRA are also shown for reference. The absorption free SED model is shown by 
a dashed line. The SEDs corrected by EBL models are shown by a dot-dashed line (LHH). 
and a solid line (LLL). 
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Fig. 3. — Electron spectrum in the emission region for the base model. The parameters for 
the injection spectrum are 7 m i n = 2, 7 max = 10 7 , p = 1.9, and the injection rate of 1CT 3 
electrons cm -3 s _1 . 
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Fig. 4. — SED of H1426+428 and a model spectrum. The model SED corrected by an EBL 
model, MLL, is shown by a solid line, and the absorption free SED model is shown by a 
dashed line. 




Fig. 5. — SED models calculated with different parameter values. Data from Whipple and 
HEGRA 1999-2000 are multiplied by exp(r 77 ) with an EBL model, LLL. Base model denotes 
the SED shown in Fig. 2 by the dashed line. 



